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We present a search for the rare decays B h^*^vT', where h^*^ stands for a hght meson. A 
data sample of 535 miUion BB pairs collected with the Belle detector at the KEKB e^e~ collider 
is used. Signal candidates are required to have an accompanying B meson fully reconstructed in a 
hadronic mode and signal-side particles consistent with a single ft'*-* meson. No significant signal is 
observed and we set upper limits on the branching fractions at 90% confidence level. The limits on 
K'^vv and B^ K'^vV decays are more stringent than the previous constraints, while the 
first searches for _B° K^vV, -k^uV, (PvV, (j)vV and — > K*'^u'u, p^w are reported. 

PACS numbers: 13.25.Hw, 14.40.Nd 



The decays B K^*^vT' proceed through the flavor- 
changing neutral-current process b —>■ svv, which is sen- 
sitive to physics beyond the Standard Model (SM). The 
dominant SM diagrams are shown in Fig. [T] Similarly, 
the decays B — *■ (tt, p)vV proceed through b — > dvV pro- 
cesses. The SM branching fractions are estimated to be 
1.3 X 10"^ and 4 x 10"^ for B K*vV and B Kvv de- 
cays [H, respectively, and are expected to be much lower 
for other modes. Theoretical calculation of the decay 
amplitudes for B — s- h^*^vv is particularly reliable, ow- 
ing to the absence of long-distance interactions that af- 
fect charged-lepton channels B — > h^*H~^l~ . New physics 
such as SUSY particles or a possible fourth generation 
could potentially contribute to the penguin loop or box 
diagram and enhance the branching fractions [1]. Refer- 
ence 0] also discusses the possibility of discovering light 
dark matter in 6 ^ s transitions with large missing mo- 
mentum. 
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FIG. 1: The quark-level diagrams for B K'vu decays. 

Experimental measurements [3| of the b ^ s transi- 
tions with two charged leptons are in good agreement 
with SM calculations [H. Further investigation of the 
forward-backward asymmetry in B — )■ K*l'^l~ Q is 



consistent with the SM although the statistics are still 
limited. Due to the challenge of cleanly detecting rare 
modes with two final-state neutrinos, only a few studies 
of K^*^v'U have been carried out to date [1,0,0]; there is 
only one examination of the corresponding b d tran- 
sitions In this paper, we report our first search for 
the decays B h^*'^vv {h^*^ stands for K+ , K^, K*", 
K*^, 7r+, 7r°, p^, and (p) using a 492 fb~^ data sam- 
ple recorded at the T{AS) resonance, corresponding to 
535 X 10^ i?-meson pairs. Charge-conjugate decays are 
implied throughout this paper. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer located at the KEKB collider fsj, and consists 
of a silicon vertex detector (SVD), a 50-layer central drift 
chamber (CDC), an array of aerogel threshold Cherenkov 
counters (ACC), a barrel- like arrangement of time-of- 
flight scintillation counters (TOF), and an electromag- 
netic calorimeter comprised of CsI(Tl) crystals (ECL) 
located inside a superconducting solenoid that provides 
a 1.5 T magnetic field. An iron flux-return located out- 
side the coil is instrumented to detect mesons and 
to identify muons (KLM). The detector is described in 
detail elsewhere ifl]. 

Candidate e+e" T(45) BB events are char- 
acterized by a fully-reconstructed tag-side B meson 
(^tag)- The remaining particles are assumed to be 
products of the signal-side B meson (Bsig). The Btag 
candidates are reconstructed in one of the following 
modes: B° ^ D^*'>-7r+, D'^*'>-at, D^*^- oi*^^ , 

B+ ^ W*^°p+, D(*)Oa+, and W*^°Di*^+. 
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The D mesons are reconstructed as Z) ^^s'^ j 

K^sTT-TT^, K^TT^TT+TT-, K+ tt' 1^- , and K+ti-tj-tt'^. The 
foUowing decay channels are included for D'^ mesons: 

KItt-tt+tt° and K- K+ ._T\ic D*-_(D*°) mesons are 
reconstructed as D'^tt^ {D'^tt'^ and D'^j). Furthermore, 
D*+ D+7, D+ K^K+ and K+K-Tr+ decays are 
reconstructed. Btag candidates are selected using the 
beam-energy constrained mass Mbc = a/ E^^am ~ P% ^^'^ 
the energy difference AE = Eb — -Ebeam, where Eb and 
Pb are the reconstructed energy and momentum of the 
-Stag candidate in the T(4S') center-of-mass (CM) frame, 
and i?boam is the beam energy in this frame. 

We require Stag candidates satisfy the requirements 
Mbc > 5.27 GeV/c^ and -80 MeV < AE < 60 MeV. If 
there are multiple i^tag candidates in an event, the can- 
didate with the smallest based on the deviations from 
the nominal values of AE, the D meson mass, and the 
mass difference between the D* and the D (for candidates 
with a D* in the final state) is chosen. We reconstruct 
7.88 X 10^ and 4.91 x 10^ charged and neutral B mesons, 
respectively. 

The particles in the event not associated with the Stag 
meson are used to reconstruct a i?sig h'^*^i'T' candidate. 
Prompt charged tracks are required to have a maximum 
distance to the interaction point (IP) of 5 cm in the beam 
direction (z), of 2 cm in the transverse plane (r-0), and 
a minimum momentum of 0.1 GeV/c in the transverse 
plane. We reconstruct if ^ {t^^) candidates from charged 
tracks having a kaon likelihood greater than 0.6 (less than 
0.4) with an efficiency of 84-91% (87-92%). The kaon 
likelihood is defined by TZk = I{C,k + C,tt), where Ck 
(£^) denotes a combined likelihood measurement from 
the ACC, the TOP, and a dE/dx from the CDC for the 
K"^ (7r='=) tracks. Pairs of oppositely charged tracks are 
used to reconstruct tt+tt" decays, with an invari- 

ant mass that is within ±15 MeV/c^ (>5 cr) from the 
nominal K'g meson mass. The tt+tt" vertex is required 
to be displaced from the IP by a minimum distance of 
0.22 cm. The direction of the pion pair momentum in the 
transverse plane must agree with the direction defined by 
the IP and the vertex displacement within 0.03 rad. For 
7^0 a minimum photon energy of 50 MeV is re- 

quired and the 77 invariant mass must be within ±16 
MeV/c^ (~2.5 ct) of the nominal 7r° mass. 

The decays Bsig — * K+vV, ir+vv, Kgi/u, and tt'^uv are 
reconstructed from single K+, 7r+, Kg, and 7r° candi- 
dates, respectively. The — » K*'^vT' candidate is recon- 
structed from a charged pion and an oppositely charged 
kaon, while B+ — > K*+vV decays are reconstructed from 
a candidate and a charged pion, or a charged kaon 
and a 7r° candidate. The reconstructed mass of the K*^ 
{K*+) candidate should be within a ±75 MeV/c^ win- 
dow around the nominal K*^ {K*+) mass. Furthermore, 
pairs of charged pions with opposite charge are used to 
form B^ p^vV candidates where the tt+tt" invariant 



mass should be within ±150 MeV/c^ from the nominal 
p° mass. For B+ p^vv, a charged pion and a tt" 
candidate are used, and a ±150 MeV/c^ mass window is 
required. A meson is formed from a K+K~ pair with 
a reconstructed mass within ±10 MeV/c^ (~2 a) from 
the nominal (/> mass. 

No additional charged tracks or 7r° candidates are al- 
lowed in the event. We select Bgig candidates using the 
variable -Eecl = £'tot — Z'rcc, where -Etot and i?roc are the 
total visible energy measured by the ECL detector and 
the measured energy of reconstructed objects including 
the Stag and the signal side h^*^ candidate, respectively. 
A minimum threshold of 50 (100, 150) MeV on the cluster 
energy is applied for the barrel (forward endcap, back- 
ward endcap) region of the ECL detector. The decays 
B D*lv are examined as control samples; the observed 
-Eecl distributions are found to be in good agreement 
with Monte Carlo (MC) simulations [13]. The signal re- 
gion is defined by -Eecl < 0.3 GeV while the sideband 
region is given by 0.45 GeV < -Eecl < 1.5 GeV. 

The dominant background source is BB decays involv- 
ing a 6 ^ c transition. A lower bound of 1.6 GeV/c on 
P* , the momentum of the /i*^*-* candidate in the Bsig rest 
frame, suppresses this background, while an upper bound 
of 2.5 GeV/c rejects the contributions from radiative two- 
body modes such as i3 — > K*^. The P* requirement is 
removed for candidates due to the lack of theoretical 
calculations for Bd ^ 4> form factors. Furthermore, the 
cosine of the angle between the missing momentum in the 
laboratory frame and the beam is required to lie between 
—0.86 and 0.95. The missing momentum is calculated us- 
ing the momenta of the reconstructed Stag and /i^*^ can- 
didates. These criteria suppress backgrounds with par- 
ticles produced along the beam pipe. Other background 
sources, such as e+e~ qq {q = u,d,c,s) continuum 
background and rare B decays involving b ^ u, b ^ s, 
OY b ^ d processes, are found to be small. 

The reconstructed -Becl distributions are shown in 
Fig. O The E^cl distributions of background are esti- 
mated with MC simulations; in particular, a large b ^ c 
MC sample corresponding to ten times the data lumi- 
nosity is introduced with a preselection on the gener- 
ator information. The background -Eecl distributions 
are normalized by the number of events in the sideband 
region. None of the signal modes has a significant sig- 
nal. Including the effects of both statistical and system- 
atic uncertainties, an extension of the Feldman-Cousins 
method 
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12l | is used to calculate the upper limits on 

the branching fractions. The observed number of events 
in the signal box and sideband region, expected back- 
ground contributions in the signal box, reconstruction 
efficiencies, and the obtained upper limits at 90% confi- 
dence level (CL) are shown in TablelH The reconstruction 
efficiencies are estimated with MC simulations using the 
B form factors from Ref. The B° (j)vV MC 

samples are generated with the B K* form factors. 
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FIG. 2: The -Eecl distributions for B h''*^vV decays. The shaded histograms show the background distributions from 
MC simulations and are normahzed to sideband data. The open histograms show the SM expected signal distributions for 
B — » K^*''uV decays multiplied by a factor of 20 for the comparison. The vertical dashed lines show the upper bound (left) of 
the signal box and the lower bound (right) of the sideband region. 



TABLE L A summary of the number of observed events in 
the signal box (A^obs) and sideband regions (Aside), expected 
background yields (Nt) in the signal box, reconstruction effi- 
ciencies including both -Stag and Bsig (e), and the upper limits 
(U.L.) on the branching fractions at 90% CL. 



Mode 


Nabs 


Aside 


Nt 


e(xlO-") 


U.L. 






7 


16 


4.2 ± 1.4 


5.1 ±0.3 


< 3.4 X 10" 


4 


K'+vV 


4 


18 


5.6 ±1.8 


5.8 ±0.7 


< 1.4 X 10" 


4 




1 


7 


2.3 ± 1.2 


2.8 ±0.3 








3 


11 


3.3 ± 1.4 


3.0 ±0.4 






K+fV 


10 


60 


20.0 ±4.0 


26.7 ±2.9 


< 1.4 X 10" 


5 




2 


8 


2.0 ±0.9 


5.0 ±0.3 


< 1.6 X 10" 


4 




33 


149 


25.9 ±3.9 


24.2 ±2.6 


< 1.7 X 10" 


4 




11 


15 


3.8 ±1.3 


12.8 ±0.8 


< 2.2 X 10" 


4 


p°vV 


21 


46 


11.5 ±2.3 


8.4 ±0.5 


< 4.4 X 10" 


4 


p^vV 


15 


66 


17.8 ±3.2 


8.5 ± 1.1 


< 1.5 X 10" 


4 


(bvV 


1 


9 


1.9 ±0.9 


9.6 ± 1.4 


< 5.8 X 10" 
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The uncertainties associated with the background 
yields are dominated by the data sideband statistics and 
MC statistics, and are summarized in Table [ill The 
possible disagreement in the i?ECL distributions between 
data and MC is checked using wrong-flavor combinatorial 



events, and an uncertainty of 0.1-2.0 events is included. 
Possible backgrounds from rare B decays are examined 
using a large MC sample corresponding to 50 times the 
data luminosity. We change the relative normalizations 
of rare B components by ±50%, and the variation in the 
background yield (0.1-1.8 events) is included as a sys- 
tematic uncertainty. 

Various sources of uncertainties are considered for the 
signal normalization and are summarized in Table IIIII 
The uncertainties in Btag reconstruction (2.0% for B'^ 
and 9.9% for B^) are estimated by comparing the yields 
of data and MC from the i?tag candidates. Systematic 
uncertainty arising from the track and tt*^ rejection is 
studied using B D'^*^(.v decays, and an error of 2.7% 
is assigned. The uncertainties in the efficiencies for de- 
tecting a or 7r° from B^ig are estimated to be 4.9% 
and 4.0%, respectively. We also vary the B h^*"^ form 
factors used in the signal MC generation according to 
the uncertainties given by the Ref. [l3|, and an uncer- 
tainty of 0.4-3.7% on the reconstruction efficiency is in- 
cluded. A larger uncertainty of 13%, which is estimated 
from the difference between the default decay model and 
a generic three-body phase-space model, is introduced 
for B^ (jjvV decays. Furthermore, the following un- 
certainties are also considered: the number of BB events 
(1.3%), tracking efficiency (1.0-2.2%), particle identifica- 
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tion (0.5-2.0%), h^*^ mass selection (0.8-2.3%), and the 
(j) — > K'^K" branching fraction (1.2%). 

In conclusion, we have performed a search for B 
vv decays with a fully reconstructed B tagging 
method on a data sample of 535 x 10^ BB pairs collected 
at the T(45) resonance with the Belle detector. No sig- 
nificant signal is observed and we set upper limits on the 
branching fractions at 90% CL. The limits obtained for 
B^ K*^vv and B+ K^vV decays are more strin- 
gent than the previous constraints from DELPHI ^5*1 and 
BaBar [T,]- The first searches for B^ K^vV, tt^vV, 
p^vv, (pvv, and B^ K*^h'h', p^uu are carried out, 
and upper limits on the branching fraction of order 10~* 
are obtained. The limit on B~^ — > tt+w is less restrictive 
than BaBar's result 0] due to a larger number of observed 
events in the signal box. The results on B ^ Ki*)i>V 
reported here are one order of magnitude above the pre- 
dictions of Buchalla et al. [l[ and hence still allow room 
for substantial non-SM contributions. A higher luminos- 
ity -B-factory experiment is required to probe the SM 
predictions for the branching fractions. 
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TABLE II: Summary of the uncertainties associated with the background yields (in the number of events). 



Uncertainties 








A"" 1/77 










(j>UV 


Sideband statistics 


1.0 


1.3 


2.6 


0.7 


2.1 


1.0 


1.7 


2.2 


0.6 


MC statistics 


0.9 


1.2 


2.6 


0.6 


2.0 


0.8 


1.3 


1.6 


0.6 


MC/data difference 


0.3 


0.3 


1.5 


0.2 


2.0 


0.3 


0.9 


1.4 


0.1 


Rare B 


0.1 


0.3 


0.2 


0.2 


1.8 


0.1 


0.1 


0.9 


0.3 


Total 


1.4 


1.8 


4.0 


0.9 


3.9 


1.3 


2.3 


3.2 


0.9 



TABLE III: Summary of the relative uncertainties for signal normalization (in %). 



Source 




K*+vV 


K+vV 












cj>vv 


N(BB) 


1.3 


1.3 


1.3 


1.3 


1.3 


1.3 


1.3 


1.3 


1.3 


Tracking efficiency 


2.1 


1.1 


1.0 




1.0 




2.2 


1.1 


2.0 


i^g/yr" reconstruction 




4.4 




4.9 




4.0 




4.0 




Sub branching fraction 


















1.2 


Particle identification 


1.3 


0.7 


0.7 




0.5 




1.0 


0.5 


2.0 


MC statistics 


3.5 


2.4 


1.9 


3.2 


2.0 


2.8 


3.3 


3.4 


2.3 


Mass selection 


0.8 


2.3 










1.1 


2.6 


2.0 


Stag reconstruction 


2.0 


9.9 


9.9 


2.0 


9.9 


2.0 


2.0 


9.9 


2.0 


Track/Tr" rejection 


2.7 


2.7 


2.7 


2.7 


2.7 


2.7 


2.7 


2.7 


2.7 


B — > /i'*' form factor 


1.6 


1.3 


2.6 


0.4 


3.0 


1.0 


1.7 


3.7 


13.0 


Total 


5.9 


11.8 


10.9 


6.9 


11.0 


6.1 


5.8 


12.5 


14.2 



